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Abstract We investigate the low-noise regime of a large class of probabilistic cellu- 
lar automata, including the North-East-Center model of Toom. They are defined as 
stochastic perturbations of cellular automata belonging to the category of monotonic 
binary tessellations and possessing a property of erosion. We prove, for a set of ini- 
tial conditions, exponential convergence of the induced processes toward an extremal 
invariant measure with a highly predominant spin value. We also show that this in- 
variant measure presents exponential decay of correlations in space and in time and is 
therefore strongly mixing. 

Keywords Probabilistic Cellular Automata • Phase Transition ■ Toom Models • 
Monotonic Binary Tessellations 



1 Introduction 

Probabilistic cellular automata (PCA) are a type of models of multicomponent systems 
whose time evolution is driven by local interactions blurred with some noise. More 
precisely, they are discrete-time stochastic processes with the Markov property, made 
of lattices of components whose individual states take values in a finite set and are 
simultaneously updated at every time step. The transition rules involve interactions 
between neighbouring components. 

The long-time limit of these processes has been the subject of many numerical 
and theoretical results in the last fifty years and a lot of questions remain open - see 
for instance the surveys [14, 15]. Despite the apparent simplicity of their discrete con- 
figuration space and merely local interactions, PCA exhibit a variety of macroscopic 

A. de Maere ■ L. Ponselet 

Universite catholique de Louvain, IRMP, Chemin du cyclotron 2, bte L7.01.03, B-1348 Louvain-la-Neuve, 
Belgium. 

Tel.: +32-10-473279 
Fax: +32-10-472414 
E-mail: lise.ponselet@uclouvain.be 



2 



Augustin de Maere, Lise Ponselet 



tendencies. Moreover, due to their basic definition, these models provide a favorable 
field of study in order to give firm foundations to the growing knowledge about non- 
equilibrium phenomena. In particular, interest concentrates on the existence of phase 
transitions in this context of non-equiUbrium statistical physics. For a non-negligible 
subset of the parameter set of certain PCA models, the limiting states of the processes 
can be non-unique and depend strongly on the initial conditions, thus providing ex- 
amples of systems which keep remembering part of the data from their remote past 
[1]. 

Here we consider the PCA resulting from small random perturbations of the de- 
terministic cellular automata in the class of monotonic binary tessellations (MBT) 
examined in [12, 13]. In an MBT, the spin variables disposed on the sites of a Z'' lat- 
tice can take only two different values and are all simultaneously updated according 
to the same monotonic function of their neighbours' states. In [12, 13], Toom gave a 
criterion for an MBT to be an eroder, i.e. to erase in a finite time any finite island 
of impurities in a predominantly homogeneous configuration of one of the two types 
of spins. The completely homogeneous space-time configuration with this spin value 
everywhere is then said to be an attractive trajectory of that MBT. He also proved in 
[13] that this erosion condition implies the stability of the homogeneous trajectory un- 
der the introduction of a small error rate: the corresponding PCA admits a stationary 
measure for which the low probability of finding an impurity is preserved. Another 
proof of this stability theorem in a more general context was later given in [3], using 
renormalization group methods. 

For example, the one-dimensional percolation PCA of Stavskaya and Toom's two- 
dimensional North-East-Center voting PCA fall within this class of perturbations of 
MBT which are eroders. They both present a phase transition as non-uniqueness of 
their invariant probability measures in the low-noise regime follows from this stabil- 
ity theorem, while, when the noise is high enough, the processes converge toward a 
unique invariant measure, regardless of the initial state [5]. 

Beside this stability result, very little is known about the long-time asymptotics of 
this class of PCA in the nearly deterministic regime. It has already been explored by 
means of simulations [1,4,9, 10, 16]. Here we address the problem from a theoretical 
point of view. By virtue of the ergodic theorem, the asymptotic behavior of such a 
PCA is actually given by its ergodic invariant probability measures. We focus on the 
statistical properties of the invariant measure with a predominant spin value which is 
obtained as a perturbation of the stable homogeneous trajectory. Reference [2] estab- 
lished exponential convergence to this stationary measure, and therefore exponential 
decay of its correlations, for the Stavskaya PCA and a class of multidimensional gen- 
eralizations of it. 

In this paper, we use a different method allowing us to extend this result to the 
North-East-Center model and all the other MBT satisfying Toom's erosion condition, 
with a restriction to MBT without memory, where the state of the system at any time 
f depends only on its previous state at time f — 1 and not directly on earlier states. 
We prove rigorously that for some well-chosen initial conditions, the process con- 
verges exponentially fast toward the invariant measure under consideration. We also 
show that this measure has exponential decay of correlations in space and in time. In 
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particular, it is ergodic and therefore extremal. The proof is based on a perturbative 
expansion, with paths and graphs, which combines a technique of decoupling in the 
pure phases previously introduced and developed for coupled map lattices [6, 8] with 
graphs constructed for a Peierls argument in the proof of stability in [13]. 

2 Formalism and Results 

We consider the cellular automata in the class of monotonic binary tessellations de- 
scribed in Section IV of [13] except that we restrict ourselves to MBT without mem- 
ory. Some of our notations are similar to the notations of [13] in order to facilitate the 
use of the graphs constructed there but they are not identical because we don't need as 
general a framework here. 

Let 5 = { — l,l}be the binary state space for a component of the system, which 
we will call a spin, and let X = S^'' be the configuration space for the whole system. 
For any spin configuration © e X, (Ox will denote the value of a at site x and (Oa will 
denote the vector {(Ox}xeA for any A C Z''. We will also use the notation (©=^,,,fl) for 
the configuration obtained from (O by replacing the spin (Ox at site x with the value 

fl e 5. 

In a tessellation, the evolution law for the state of the spin at a site x of the Z'' 
lattice involves its neighbours, which are defined as the elements of U{x) = x-\-U 
for a fixed finite set [/ = {mi , . . . , C Z''. Starting from a configuration (O ^X, the 
simultaneous updating of every spin at every time step consists in transforming the 
spin (Ox at site x into the value 

^x{(0) = <j){C0x+ui,---,C0x+UR) 

for a given function (j) : S'^ ^ S. From now on, as we consider monotonic binary 
tessellations, (j) will be monotonic in the sense that if a)„ < (0^ for all u G U, then 
0(®f/) ^ 0(^c/)- We also reject the trivial case of a constant function (p. Note that 
these two assumptions imply that 0(±1,...,±1) = ±1 and then that the configurations 
and { — l},£2rf are left invariant by the deterministic time evolution. They 
generate completely homogeneous trajectories. 

In addition, we suppose that (j) verifies the Toom erosion criterion. It is expressed 
in terms of the plus sets: the subsets Z C [/ such that if Ct)„ = +1 for all u G Z, then 
0(^t/) = +1- When applied to the particular case of an MBT without memory, the 
criterion given in Proposition 1 of [12] and in Theorem 6 of [13] states that the inter- 
section of all convex hulls of plus sets in R'' must be empty. 

We now describe a way to introduce some noise in such an eroding MBT to gener- 
ate a PCA. The system follows the same updating rule as in the deterministic case but, 
at each site of the lattice and at each time step, an error can occur with a probability 
of order e, for some e in [0, 1], the spin then taking the opposite value. Occurrence of 
an error at a site is independent from occurrence of errors at other sites. The process 
resulting from the sequence of simultaneous updatings of all spins becomes a stochas- 
tic process. The local transition probabilities p{£,x\^) of that process only depend on 
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^u{x) ■ We suppose that they satisfy the following assumptions, for some given values 
of the parameters e G [0, 1] and a > 0: 

(Al) if ^ then p{^^\a) < e; 

(A2) if fl = then \p{^x\^) - p{ix\^^y,a) \ < ap{^x\m) for ally e Z''. 

Assumption (Al) is the low-noise condition: it ensures that, at each site of the lattice, 
the deterministic rule is followed with a probability of at least 1 — e. Note that we 
make no restriction about the possible bias in favor of errors of a specific spin sign. 
Assumption (A2), used in conjunction with the monotonicity of (j), ensures a form of 
decoupling in the pure phases: if a spin was flipped without changing the deterministic 
prescription given by the value of 0^, it would be of little consequence to the involved 
transition probabilities. The local transition probabilities define a transition probability 
kernel on X, which can be formally written as 

p[ii®] = n /'(^-^i®)- (1) 

Assume that S is endowed with the discrete topology and X with the product 
topology. Let "^(X) be the set of continuous functions fromX into IR with the norm: 
l^l^ = supi^^^ \(p{^)\. Since X is compact, "^(X) with the norm |.|^ is a Banach 
space, and its dual is (X), the space of signed Borel measures on X with the norm: 

\^\ = sup{^i(p)\(pe'S'{x),\(p\^<i}. 

For PCA on infinite spaces, it is often useful to consider the following semi-norm on 
'^{X): if 5x(p{a) = (p{a) - (p{a^x, -(Ox), we define: 

ll'Pll= L |5v(PL. (2) 

xezd 

Let T : ^{X) — > {X) be the transfer operator associated to the transition prob- 
ability kernel of (1), explicitly defined as: 

Tfi{(p) = |d«|d|9(|)P[||«]M(«), (3) 

for any (p m'f{X). 

Let 1_ A denote the indicator function of the subset {a E X\cOx = — 1 Vx G A} 
and, for any measurable set Y CX, let ly ■ ^^{X) — s- ./^ {X) be the operator defined 
as: 

=^((plr) V(pG'if(X). 

As the MET defined above satisfies the Toom erosion condition, the results of 
[13] apply to it. Theorems 5 and 6 of [13] imply the stability of the homogeneous 
trajectory with the spin value +1 everywhere. More precisely, let 5'^' be the Dirac 
measure concentrated on the configuration {+1 l^g^rf- Then, using these theorems and 
Definition 1 of a stable trajectory in [13], we have 

lim sup r"5<+' (1_ f.i) = 0. (4) 

xez'' 
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Now consider the following sequence of measures: 




We can always extract from it a weakly-* convergent subsequence. Indeed, the Banach- 
Alaoglu theorem states that the unit ball of ^(X) is compact in the weak-* topology. 
The associated limit is explicitly given by 

= lim [- E r*5W((p)) V(pe<r(X), (5) 



for a certain subsequence [rij) of increasing positive integers. The measure yi^^^ is 

obviously an invariai 
tion (4) implies that 



obviously an invariant probability measure, in the sense that T = [l^^^, and equa 



lim sup Ai^- (!_,{,}) =0. (6) 

Example 1 (The Stavskaya model) The percolation PCA of Stavskaya [11] is one of 
the first PCA for which the existence of a phase transition has been rigorously proven 

- see an outline of the proof using the method of contours in Chapter 1 of [15]. It has 
c/ = 1, f/ = {0, 1} C Z and 0((Oo,ft)i) = -1 if and only if coo = COi = -1. Errors are 
strongly biased in the sense that they can only transform a spin +1 into a —1 but not a 

— 1 into a +1. The Dirac measure 5'"' is then always a stationary state. We will call it 
j'^inv- Besides, the Toom criterion is satisfied by these U and since the intersection of 
the plus sets {0} and {1} is empty, therefore equation (6) carries the existence of an 
Ec > such that for all e < Ec, the process admits a second different invariant measure 

^^lnv 

Example 2 (The Toom model) The North-East-Center voting model has d = 2 and 
U = {(0,0); (1,0); (0, 1)} C 1? and returns the majority spin among these three 
neighbours' states. This model verifies the Toom erosion condition too, because 
{(0,0); (1,0)}, {(1,0); (0,1)} and {(0, 1); (0,0)} are plus sets whose convex hulls 
have an empty intersection. Moreover, the model is symmetric under the interchange 
of positive and negative spin signs. Combining this symmetry with definition (5) gives 
two invariant measures, jx^^^ which satisfies equation (6) and /Xj'jj^ which satisfies the 
symmetric counterpart of equation (6), thus revealing that they differ as long as e is 
small enough. 

Then, in such models, if e is below a critical threshold e^ , there exists an infi- 
nite number of invariant probability measures, since any convex combination of pL^^^ 
and /ijjj^ is also an invariant probability measure. But, among these invariant mea- 
sures, some have interesting properties such as ergodicity and exponential decay of 
correlations. In [2], the low-noise regime of a class of PCA is examined, including the 
Stavskaya model and a multidimensional generalization of it, but not the Toom model. 
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Exponential convergence toward equilibrium of the processes with the initial condi- 
tion 5'+' is proven for these PCA by constructing a cluster expansion. Exponential 
decay of correlations in space and in time follows for the non-trivial natural invariant 
measure jx^^l that appears at the phase transition of the Stavskaya model. 

Our argument extends the results of [2] to the whole class of PCA defined above. 
It relies on an expansion which isolates the influence of each space-time site on each 
other site in its future, along several paths of influence. These paths pass through 
the dominant phase with positive spins almost everywhere and Assumption (A2) of 
decoupling in the pure phases can be used in order to bound the influence of spins one 
on another. The paths seldom encounter a negative spin. Whenever they do, in view to 
evaluate how improbable that negative spin is, we will attach to it a graph borrowed 
from a Peierls argument in the proof of stability in [13] and bring Assumption (Al) 
into play. The fact that paths select one site at a time will allow us to make use of the 
one-dimensional Toom graphs by associating them to single sites separately. As for 
the cluster expansion in [2], it requires contours which can enclose blocks of several 
spins, as the Stavskaya contours do, but not the most general Toom graphs, as far as 
we can see. 

In Sections 3 to 6, we will show that any initial probability measure in a suit- 
able basin of attraction of ^(X) converges exponentially fast toward jl^^l- We will 
consider the sets 



with K >0 and e' G [0, 1] and prove the following result: 

Theorem 1 For any MBT without memory characterized by a non-constant mono- 
tonic function <j) verifying the Toom erosion condition, there exist a* > and a posi- 
tive and strictly decreasing function defined on [0, OC* [ such that, for all a £ [0, OC* [ 
and £ G [0, e*(o£)[, the following assertion is true for any PCA satisfying the corre- 
sponding Assumptions (Al) and (A2). For any probability measure /X in ^^^\K,e') 
with K >0 and e' < e* (o:), there exist some constants C <°° and C7 < 1 such that, for 
all (p e'^{X) and all n e N, 



Note that for PCA which are symmetric under the interchange of the two types of 
spins, such as the North-East-Center model, the symmetric result for in the class 
{K,£') and for ^^^l is also valid. 
Eventually, in Section 7, we will prove that ji^'^l has exponential decay of correla- 
tions in space and in time and is, consequently, strongly mixing. 

Corollary Assume that a < and that e < e*(o;), with a* and e* as given by The- 
orem 1. Then there exist some constants C' <°° and r\ < I such that for any (p, \j/ in 
^{X) with \\(p\\ <oa and || <°°, and with a positive Manhattan distance d{(p,\^) 
between their supports, we have 
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For the exponential decay of correlations in time, we define the operator T : '^(o{X) 

r<p(o))= /d^(p(^)p[^|«], 



which is simply the dual of the transfer operator acting on signed measures. 

Theorem 2 Assume that a < OC* and that £ < £*(«) with a* and e* as given by The- 
orem 1. Then, for any continuous functions (p, \j/ with finite supports, there exist some 
constants C(p^\j/ < °° and (7 < 1 such that, for all n 6 N, 



Mi^!('pr'v/)-M"((p)M"(v/) 



3 Path Expansion 

In this section, we will introduce a path expansion which is essentially equivalent to 
the Dobrushin criterion [5], using here a formalism which was originally introduced 
by Keller and Liverani for Coupled Map Lattices [6]. 

Let -< be any well-ordering of Z''. The operator iTj : '■^ {X) — > '■^ {X) is defined as: 

where, from now on, a will denote the configuration for which a^ = a = +1 for all 
X £ Z"^. {^yx,a^x) is the configuration obtained from a by replacing the spins at 
all sites y ^ x with the value a. Using telescopic sums, we can check that, for any 
continuous function (p, 

(p{m}^(p{q)+ nM^)- (7) 

xez'' 

With a slight abuse of notation, let us denote by Fix : ^{X) ^{X) the dual of 
the operator Tlx '■ "^O^) '^(^)- The image of ^{X) under this operator is actually 
included in the set 

^x = {/i G {X) I ^(9) = if (p{(0) is independent of (Ox) (8) 

whose elements verify the following property: 

jue^, |ju((p)| < |ai||5v(pU. (9) 

With equation (7), we can see that any signed measure of zero mass pL S ^{X) 
with =0 admits the following decomposition: 

Ai= L ^-Ai- (10) 

A-eZ'' 

While Tlxii belongs to ^x, it is no longer the case for TTIxfi. Nevertheless, since the 
interactions are local, we will see that TTIxfi can be expressed as the sum of R signed 
measures: a first one in ^x-u^, a second one in ^x-u2^ ■■■ ^ and a last one in ^x-ur- 



Augustin de Maere, Lise Ponselet 



For this, consider an arbitrary measure jix G Using definitions (3) and (8) of T 
and together with our hypothesis that p{£,y\w) only depends on ffls+j/, it is easy 
to check that 

rAi. = ^r(--"'-'^Vx, (11) 
(=1 

where R new operators have been defined: 

r(^-"""V(<P)=/do)y"d|(p(|)T("-""^)(|,«)Al(«) fori=l,...,R (12) 
with the kernels 



We notice that, for all /, the image of .^^(X) under the operator T^^ ""'^^ is included 
in ^x-ur 

Let now jihea signed measure of zero mass and consider r"/i. Using the decom- 
position (10) and applying (11) iteratively, we find 

L L ••• E T^"-''"-'^---T^''''°^nx,^. 

This sum can be rewritten as a sum over paths. Indeed, if we introduce 
Px=[r:{0,...,n}^r' I 7,=xand7,_i-7, GC/Vf}, 
it is equivalent to the following compact expression: 

T"^=L L rir'^'^^'-''iT,„M, (14) 

where the operators T^^''^'-''> have to be applied in chronological order 

In the case of a weakly interacting system, this sum over paths can be used to prove 
the existence of a unique invariant probability measure, under the assumptions of the 
Dobrushin criterion [5]. The system we are considering here is certainly not weakly 
interacting. However, in order to prove Theorem 1, the idea will be to take advantage 
of the decoupling in the pure phases conveyed in Assumption (A2). This assumption 
will provide upper bounds of order a on the couplings r*^'-^' '' in the positive phase. 
Indeed, as a = +1 in (13), those bounds will be obtained for the instants t such that the 
considered space-time configuration presents at time r — 1 the value +1 for 
As for the negative phase, it will be shown to be infrequent enough, for a suitable 
choice of initial condition. 



1-1 



{p{ix-u,\o)) - pii,-,.\a^x,a)^ f n Pi^-^-"j\^)] ■ (13) 
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4 Pure Phase Expansion 

For fixed x E Z'' and 7 e P^, if 7 denotes both the function and its trajectory 
{{t,Yt) 1 1 =0,...,n}, let us partition the trajectory into 7'+' C 7\ {70} and 7*^ = 
(7\ {70}) \ 7<+'. We define, for f = 0, . . . the sets F C X with the follow- 

ing indicator functions: 

If n ^+iuo})) n ^-iuo^))- 

.v;(r+l,.v)G7(+) .r:(^+l,;c)Gy(^ 

These subsets lead to a partition of the configuration space X^" ■ ■ "^ Inserting this 
partition in ( 14), 

T"^^=L L I ri(V(rW.o^*'''''-'')W'.o)^roM, 

where 1^ is nothing but the identity operator. Using property (9) together with 
the fact that the image of ^{X) under r^^ ''" ') is included in we have, for all 

\T"^m< ILL 

xeZd y£-Pv7(+)Cy\{yo} 

(15) 

Next, for given x e Z'', 7 e and 7*^' C 7\ {70}, we define the set 
<f (7*^') = {{bYU e X^'-"^ \B^F (7<+>,f) Vf e {0, ...,«}} . 

We will now build a graph using the construction in [13] in order to control the extent 
of the negative phase. We will then rewrite the expansion (15) as a cluster expansion 
in terms of a combination of paths and graphs. 

5 Graphs 

As we included the Toom erosion condition in the hypotheses of Theorem 1, we know 
that the intersection of all convex hulls of plus sets of the considered MBT is empty. 
Using Lemma 12 of [ 13], it implies the existence of an integer q 6 No, of q linear func- 
tions Li,... ,Lij from IR''+' to IR and of a positive r such that the hypotheses of Theorem 
1 of [13] are verified by the MBT and by the homogeneous space-time configuration 
with positive spins everywhere. In particular, we can use the structures which are de- 
fined in the Peierls argument proving this stability theorem in Section II of [13]. We 
will not rewrite the whole construction here but only extract the objects of use for our 
purpose and their properties. We will also adapt the general formalism to the more 
restricted class of models under consideration here. We refer the reader to [13] for 
all details and, for an introduction to that very general proof, to the review given in 
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Appendix A of [7] in the particular case of the North-East-Center model. A good ex- 
ercise could also be to consider the case of the Stavskaya model. One then obtains the 
contours described in Chapter 1 of [15]. 

In order to bound the probability of finding a negative spin at a given time f_ and 
at a given site x_ of the lattice, when the initial measure is 5*^', [13] associates a 
structure called truss to each space-time configuration presenting this negative spin. 
Its iterative construction is given in Section II, Part 3 of [13] where it is denoted 
by n. This truss consists of an unordered sequence of edgers, which are maps from 
{ 1 , . . . , ^} to {0, . . . , f_ } X J}' whose image contains exactly two points. Here we use, 
much like in the review of [7], an equivalent reformulation. We view an edger as an 
unoriented edge that links two points of {0, . . . ,f_} x Z'' and bears an extra attribute, 
namely a partition of the set {1, . . . of poles between its two vertices. This natural 
bijection between edgers and edges extends to trusses which we treat as graphs on 
{0, . . . ,f_} X Z'' composed of these edges corresponding to their edgers. We name 
them Toom graphs. 

The Toom graph G thus associated to the space-time configuration identifies part of 
the errors which lead to the negative spin at (f_ its set of vertices Vq admits only 
points where the given configuration presents negative spins and it contains (f_ , x_ ) . It 
traces the source of this negative spin, via their propagation according to the determin- 
istic rule conveyed by the function 0, back to some previous errors (when 0.t(c^-i) 
= +1 but (£)t,x = —1)- Such errors, classified in a distinguished subset Vg of Vc, all 
carry a probability smaller than e, because of Assumption (Al). Vq is noted S in [13] 
and constructed iteratively parallel to the truss. The large number of different graphs 
associated to all possible space-time configurations with a negative spin at (f_,x_) is 
then an 'entropic' factor which has to be counterbalanced by an 'energetic' factor, i.e. 
the probability of the graphs, which is bounded by the probability of making all the 
identified errors, g'^ol. 

Now, using our reformulation of trusses in terms of Toom graphs, the construction 
algorithm of the trusses in Section II, Part 3 of [13] guarantees the following proper- 
ties: 

(PI) The graphs are connected and contain (f_,x_), which is their only ver- 
tex at time coordinate f_. We will call it their origin. Moreover, their 
edges are either arrows, with a displacement vector between their ver- 
tices of the form ±(— 1,m) for a m in f/, ot forks, with a displacement 
vector (0,m''' — m*^^^) where m'^^ and m'^^ are both elements of U. There- 
fore, taking into account the 2"? possible allocations of the q poles to the 
two vertices of an edge, for any given point of {0, . . . ,f_} x Z'' at most 
21 {R^ +2R) different types of edges can have that point as vertex. For 
any such connected graph, there always exists a walk which starts from 
(f_,x_), passes along every edge exactly twice and then comes back to 
its departure point. We choose such a walk and consider the sequence 
that records, at each of its steps, the displacement vector and the pole 
distribution of the travelled edge, and call it the Eulerian walk associated 
to the graph. This correspondence is injective. Consequently, the num- 
ber of different graphs grows only exponentially with their number of 
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edges, \Eq\, as the number of Eulerian walks of length 2\Eq\ is less than 
[2''(/?2 + 2/?)]2|£ol. This is Lemma 3 of [13]. 
(P2) The number \Eg\ of edges is in turn related to the number \Va\ of identi- 
fied errors: 

YT^\^o\ + ^<\Va\- (16) 

This is due to a principle which presides over the graph construction in 
[13] so as to take advantage of the erosion property: the truss is overall 
even and the number of forks equals the number of identified errors | Vq \ 
minus one. Using Lemma 2 of [13], one deduces the upper bound (16). 

Reference [13] uses these properties to prove that the 'energetic' factor wins as long 
as e is small enough. The probability of finding a negative spin at a given time and a 
given site can be made as small as desired by taking £ appropriately small. 

Here in our treatment of the space-time configurations that present, for all t such 
that (f , 7r ) G 7*"', a negative value of the deterministic prescription given by evaluated 
at time f — 1 in 7, + [/, we will hardly need to modify Toom's construction, but we will 
associate to each space-time configuration a collection of Toom graphs instead of only 
one, as we are interested in the source of a collection of negative spins, corresponding 
to the different elements of 7* ' . 

We define a map 

g:S{t'')-^'^{t')=g{'S{t'')) (17) 

where G is a graph made of a disconnected collection of Toom graphs whose con- 
struction, for a given (ffli)JLo S (7*^'), consists in the following steps. 

1. We pick fi, the largest t such that {t\ , 7,1 ) G 7*"'. We know that ^y,^ (o)ri-i) = -1- 
We consider (w,)"^o which is obtained from by replacing Oiti,y^ with the 

value — 1 . We construct the Toom graph G\ having the negative spin at (f 1 , 7| ) as 
origin and which is associated to (w,)JLo according to the algorithm in Section 
II, Part 3 of [13]. The whole algorithm remains valid here, the only difference be- 
ing that no assumption about the initial condition forbids the presence of negative 
spins at f = in ©0 while, in [13], the initial condition is 5'+' and therefore the 
graph is always contained in { 1 , . . . , f_ } x Z'' since the points where the configu- 
ration has positive spins can't be vertices. A rule is then added to the algorithm: 
when the graph under construction reaches negative spins at f = 0, they are con- 
sidered equivalent to errors and classed as elements of Vgi ■ In other words, this 
amounts only to a one-unit shift of the initial condition 5'+' along the time axis. 
So we end up with a connected graph G\ . Its set of vertices VQ^ contains {t\ , y,^ ) 
and at all points (f,jc) G Vci \ {(^b 7ri )}> ths space-time configuration (oj()"^o takes 
the value (Ht.x =■ — 1 . It has a distinguished subset of vertices Vci Q ^Gi where the 
space-time configuration presents errors: 

V(f ,;<:) G Vbi such that t ^ 0, = +1. 
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It has a set Eq^ of edges connecting its vertices and satisfying (16). 
Gi is the first of the Toom graphs whose union will form the graph G. 

2. We pick the next maximum f2 < f i such that (f2 , ) € 7*"' . We perform exactly the 
same construction as before and associate to (f2, /rj) the graph G2 with properties 
analogous to those of G\. If Vg, H Vg, 7^ then we discard G2, as we can't count 
any error twice. Otherwise the union of Gi and G2 will be part of G. 

3. We repeat the same process up to the lowest time t such that (f , 7r) £ 7* discard- 
ing any Toom graph G,- which intersects one of the previous retained ones. 

4. G is defined as the union of all the retained G,, / G { 1 , . . . , 1 7* ' | }. is the union 
of the retained Vg, , Vg is the union of the retained Vg, and Eq is the union of the 
retained Eq. . 

The map ^ of (17) is then completely defined. It is usually not injective, as there are 
many configurations corresponding to one graph. But S" (7*^') can always be written 
as; 

S{t')= U 8''G. (18) 
Gef/(7(+)) 

Defining, for all G e ^ (7*+)) and for f = 0, . . . ,n, the subsets E{G,t) C X: 

l£(G,r)(®)= n n (19) 

a: (r ,A-) e l^oA r'"' : (' + 1 ,-v) e Vg 

we have, by construction of the map g, 

V'G((®0;U)<rili?(G.r)(^). (20) 

r=0 

Let us consider any G G (7*^') • G is the union of c individually connected but 
pairwise disconnected Toom graphs (G,| , . . . ,G,-^,). From this point on, they will be 
noted (Gi , . . . , G,-, . . . , G^). They all possess Properties (PI) and (P2). We show that G 
itself inherits similar properties. 

First we find an upper bound on the number of different graphs with given numbers 
of edges and of connected parts. We define a map w on (7*^'). For G G ^ (t*^')' 
w(G) consists of two elements; first, the list of the origins of Gi, . . . ,Gc; second, the 
sequence given by the concatenation of the Eulerian walks associated to Gi Gc in 
Property (PI). This map w ; (7*^') — > w (^^ (T*^')) injective. 

Proof The induction procedure described in Section II, Part 3 of [13] for the construc- 
tion of trusses starts from the vertex which we called the origin of the graph. One can 
check by inspection of the induction steps that the first step creates exactly q arrows 
leaving from this origin and then no other edge with this vertex will be drawn during 
the following steps of the construction. Consequently, we know that the Eulerian walk 
associated to a Toom graph G, will leave from and come back to the origin of G, ex- 
actly q times. Given any element of the set w (7'^''))' unique inverse image can 
then be deduced from the list of origins and the sequence of steps of the concatenated 
Eulerian walks by the following method. Starting from the first origin recorded in the 
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list, we add the steps of the sequence and redraw the first connected part of the graph, 
until the origin has been reached q times. Then we jump to the next origin recorded 
in the list and read on the sequence of steps until we again come back q times to 
this second departure point, and so on, until we have read the whole sequence. While 
redrawing a graph, we give back each travelled edge its pole distribution, which is 
recorded in the sequence. 

For all c e IM and all m e N, let (7*+',c,m) be the subset of (7*^') consisting of the 
graphs with c connected parts and m edges exactly. Any element of w (7*+',c,m)) 
is made of a list with c origins chosen among the elements of / ' and of a sequence 
with exactly 2m steps. Therefore Property (PI) is transferred from Toom graphs to 
graphs G G (t*^' , c, m) : 

|^(7<+',c,m)| < (^'^^^\y[2''{R^+2R)f" <2\'>^'^\.[2i{R^+2R)f". (21) 

On the other hand, every connected part verifies Property (P2), that is to say, satis- 
fies the analog of (16). Summing this inequality over all parts gives a similar property 
for G: 

YTh^\^o\+c<\Vo\. (22) 

We also have 

ly-'l <|£gI+c. (23) 

Proof It is sufficient to construct an injective map / : 7*"' — ;> Eq U {1, . . . ,c}. Let us 
consider any {t,Yt) S 7*"'- If (f , 7r) belongs to a connected part G, of G which is con- 
tained in {0, ... ,f} X Z'', then we know that (f , 7r) is the unique origin of G;. We as- 
sign to (f , 7() the image /((?, 7)) = / G {1, . . . ,c}. Otherwise, we know that the Toom 
graph with origin (f,7) has been discarded. Then that Toom graph intersects at least 
one of the retained Toom graphs with origin (^,7), s > t. Now the discarded graph 
is contained in {0, . . . ,f } x Z'' and the time component of edges of Toom graphs has 
maximum absolute value 1, as can be seen in the definition of arrows and forks in 
Property (PI). Therefore there exists an edge of the retained graph which arrives onto 
some point (f x G Z'' . Such an edge belongs to Eo and we take it to be the image 
of (f , 7) under /. The map / thus defined is easily seen to be injective. 



6 Exponential Convergence to Equilibrium 

We now combine the collections of Toom graphs introduced in Section 5 with the 
paths of influence of Section 3 and their pure phase partition of Section 4. Inserting 
partition (18) in expansion (15) yields a sum over graphs. It introduces intricate cou- 
pUngs between the configurations at different times, due to the complexity of the Toom 
graph construction algorithm in Section II, Part 3 of [13]. Since 1^,-1(3 does not have 
the property of factorization over time, we use the upper bound (20) in the equivalent 
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form ((oj,)"^q) = Ig-iG ((cJi)"=o) -11^=0 l£(G.r)(^)' indicator functions can 
only take values or 1. Now, for all (Wf)JLi' we have 



r(ri-ro; 



F(7(+),0)n£(G,0)V 



(•) 



where |ny(,/i| is not the norm, but the absolute value of the measure Tly^ji. Conse- 
quently, if we define the operator f'''^^ by replacing t'^^''^^(i^,oj) with its absolute 



value 



in (12), the graph expansion gives 



(24) 



< I 

Gef#(y(+)) 



n 1 



(=1 



'-F(y(+),On£(G,0 



fiYi.rt- 



F {■/.+') fl)nE{Gfi) r^ro 



This expansion will be the starting point of the proof of Theorem 1 . Before, we 
need the following lemmas. 

Lemma 1 For all x e Z'', y E Px, 7*^' C 7\ {70} and G G ^ (7*^')' V £ 

and for all t ^ {I,. . . ,«}, 



^F{y(+lt)nE{G.t)^ ^F{y(+\t-l)nE{G,t-iy 



< 



^F(yi+lt-l)nE{G.t-l)'^ 



where we introduced the notation 7'^' = 7*^' n { (f , jc) |x e Z'' } and the analogs 7' ' and 
Vg,,- 

Proof Using definitions (12) and (19) of r'^'^^' i' and l£(G,r)j we akeady have 

n , > f(rf. rf-i)-ii ,,. V 

^F(7(+)/-l)n£(G,r-l)^ ■ 



< 



(25) 



sup 

fflGf (y(+)j-l)n£(G,r-l) 



d| n 1- (^-^ 

x:(?,x)GyG\r'^ 



r(rf,y<-i) 



Now, keeping definition (13) of t(''''1''-i) in mind, the contributions of the spins t,^ at 
different sites x S to the integral in (25) are decoupled and factorize. We first con- 
sider the set {x £ Z'^ : {t,x) €Vq \ 7}. Because of Assumption (Al), its elements all 
contribute by a factor bounded by £, since the supremum in (25) is taken over config- 
urations (O in E{G,t — 1). Everywhere else on Z'' \ {7}, the contribution is trivially 
bounded by 1. For x = 7, we need upper bounds on \p (^y, — p (^y, ,«) | ■ 
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- if (py,{a) = +1, 

- and if ^y, ^-l,\p {^y, \a) - p {^y, |«^y,_, ,«) \ < a p {t,y, \(o) < a e; 

- and if ^y, = +1, \p {^y \a) - p {^y, \(o^y_,,a) \<ap {^y, \(o) , 
where we used Assumptions (Al) and (A2); 

- if (j)y (a) = — 1, we use the trivial bound \p (^y, — p {^y |^^y,_, ja) | ^ 1- 

Therefore we obtain the following three types of upper bounds on the contribution of 
the spin at 7, to the integral in (25): the contribution is bounded by a e if (f,7r) £ 
7<+' n Vg, by a if (f, 7,) g 7*+' \ Vq or by 2 if (f, 7;) e 7*"' . Inserting all these bounds in 
(25) proves Lemma 1 . 

Lemma 2 For any MBT without memory characterized by a non-constant monotonic 
function (j) verifying the Toom erosion condition, there exist a* > and a positive 
and strictly decreasing function defined on [0, OC* [ such that, for all a G [0, OC* [ and 
£ G [0, e* (a) [, the following result is verified for any PCA satisfying the corresponding 
Assumptions (Al) and (A2). Let K > and e' G [0,e*(a)[. Let jl G ^{X) be such that 
its absolute value is in the class ^'^'(/TjC'). Then, there exist some constants 
C <°° and (7 < 1 such that, for all (p G "^^(X) and all n G N, 



LIE 

xeZd y^Px yM Cy\{yo} 



ri(V(r(+),0^'"^''-'')lF(yM,o)-f7roM 

r=l 

Proof Multiple uses of Lemma 1 on the RHS of (24) imply 



|5v(pL<C||(p||a". 



n 1 



1=1 



^F(7(+),r)n£(G.f) 



'-(Yi-ri-l 



F {■/+'> fi)nE{G,0) l^^ro 



<eE;L,I^G,,l alr'+'l2lr<- 



^F(y(+),0)nE(G,0) l-^yo^l 



The important point is now that we assumed that \ belongs to ^'+'(/r,£'). Indeed, 
since 1e{g.o) < '^-,Vg.o l-,i/oo(^tro'«^ro) ^ l-,0oo(^)' and using the definition 
of Hyg, we obtain 

V(7(+).0)n£(G,0) l-^ro^l ^ 2Ke'\^°-°\. 
Consequently, with e — max{e,e'}. 



n 1 



t=i 



(y(+',r)n£(G,r 



fiYi.rt- 



F(7(+),0)n£(G,0) 



<2K£ 



Vg\ 



a 



yW|2|yH| 



(26) 



Now for all G G ^ (t*^') with c connected parts, (23) implies that the number of 
edges can be written as \Eq\ = | y^^' \—c + k with a certain A: G N. Therefore ^ (7*^') = 
UcGN (7*^* , c, 1 7*"' I — c + fc) and, by virtue of the graph properties established above. 



ken 
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we have, using first equation (22) and then equation (21), 

G6^(y(+)) f=0^=0 



< £ £2l''"l.[2'?(7?2 + 27;)]2(lr'-'l-+^)eTTi-T(l^"l+2^'-"''-+^) 

c=0 k=0 

Ir'-'l 

2.[2''{R^ + 2R)]^e'+^-"- 

< 

2c,,-- 



1 - [2'i{R^+2R)Y~e 1+2'"-' 1 - [2'=/(/;2 +2/;)]-2£ i+v 



provided that £ and e' are such that [2'i [R- + 2R)]'^e^l^^+^'i'' '> < 1. 
Combining (24), (26) and (27), we find 



7^ 



(27) 



LEE 

j:GZ''ye^v7(+)cy\{yo} 



r=l 



1 



<CEE E (4.[2'?(7;2+27;)]2£i+2.'-' ) \5,<p\ 

x-eZ''ye/'xy(+)Cy\{yo} 



where 

2/: 



C = 



1 - [2'?(7;2+2/;)]2£i+2,,-i 1 _ [2'?(7;2+27;)]-2£i+2<?'' 



Here, we can use Newton's binomial formula: for any finite set A and any x and y in R, 
Ibca-^'^'j''*^^' = {x+y)\'^\. Finally, since = R" and keeping definition (2) in mind, 
we get the desired upper bound if we take CJ = R{a + A.[2''{R^ + 2/?)]2£ V(i+2<?'-"')), 
For a to be lower than 1, the parameters must satisfy a < a* = 1//? and £ = 
max{£, £'} < £* (a) where the function £, is strictly decreasing and positive on [0, a, [. 

Using Lemmas 1 and 2 and the previous sections, we now prove Theorem 1 . 

Proof (Theorem 1) We take the constant a* and the function £, obtained in Lemma 
2. Let us consider any PC A satisfying Assumptions (Al) and (A2) with parameters 
a < a* and £ < £*(«). We write T"pi{(p) - Mini(^) = T"p.{(p) where il= fx- n'^l e 
^{X) is a signed measure of zero mass. Therefore we can apply the above path 
expansion and pure phase expansion to T"ix which then satisfies (15). Using defini- 
tion (5), pi can be rewritten as /i lim*^„ ^^^i^ (jU - so that \T"pi{^)\ is 
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bounded above by 



I "7-1 



LEE E 



t=l ^ ' ' 



1 k=0 xeZd r^Pi y(+) C 7\ { yo } 



(28) 



As the probability measure jx belongs to iK,£') with e' < e*(a), Lemma 2 
applies to the first part of (28): 



E E E 



t=l ^ ' 



|5,(pL<C||(p|1cj«, (29) 



where C <°° and C7 < 1 are given in Lemma 2. 

As for the second part of (28), it can be bounded thanks to a slight modification of 
the same arguments. We don't know whether all considered models have the property 
that r^^S'^' belongs to some {K,£'). Therefore, we will now extend the Toom 
graphs up to time —k instead of 0. For fixed A; G N, x e Z'', 7 e and 7*^' C 7\ {70}, 
we define 



{-k n} 



g),eF{y-+\t) Vf e{0,...,«}}. 



A map gk is defined, similarly to the map g above: 

(.m)"=-k-^G^gk{im)"=-k)- 

The only change in the graph construction algorithm takes place whenever a Toom 
graph Gi under construction reaches a negative spin at time t = 0. Instead of classing 
it into the set of identified errors Vg, - we carry on the construction of this branch 
of Gi, as for positive times, until we meet either an error or a negative spin at time 
—k, which is now considered the initial time, and class it in Vq.. Again, this merely 
amounts to a translation of the initial condition 5'^' by k units along the time axis. 
At the end, we obtain a disconnected union of Toom graphs on {—k, ,...,«} x Z'', 
with origins (f , 7) G 7* ' and with the properties described above. In particular, all the 
previous results about (7*^') still hold for (7*^') and for the associated subsets 
^/t (7*^' 1 c, m) of graphs with c connected parts and m edges. 

Extending definition (19) of E{G,t) to graphs G G J^/t (7*^') and to negative times, 
the analogs for gi^ of {\ 8) and (20) lead to the following graph expansion for the second 
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part of (28): 



(+) 



< E 

GG'4(r(+)) 



n 1 



'-F(7<+),0n£(G,r) 



(30) 



°lF(y(+),0)n£(G.O)®ro n (l£(G,0^)li?(G,-A')5<^' 
r=-*:+l 

where we defined the operator 0^: 

0,M(<P) = /d« /d^ (p(^) (iT,P[.|«]) f<^) 



with JTv acting on the measure P[- 

Lemma 1 as well extends immediately to graphs G G (t*^')- ^PPlyii^g i'^ " times 
to the RHS of (30), we find that it is bounded above by 



^ £i;L,l%.l«lr«l2lr"l 

GG;ft(y(+)) 



V(y(+).0)n£(G.O)®ro n (l£(G,r)r)l£(G,-A')5'^' 



Then we can use again Assumption (Al) k times. Indeed, we know from definition 



l-yG/l£(G,r-l)V 



< 



(19) that l£(G/) !i -"--Vgi from Assumption (Al) that 

gl^'o.rl |l£((j,_i)v| for any probability measure v. The operator 0^ can be handled 
in the same way as T, taking its definition into account, together with the fact that 
^-Vqo i^hx J -<.r ) !i 1 _ Vq Q (.^ ) ■ It will simply introduce an extra factor of 2 due to the 
operator iTyg. Lastly, by definition of 5'+^ 5'^' {^E{G-k)) = unless Vc-k is empty, 
in which case 5'^' (l£:(G,-/t)) = 1- Consequently, the RHS of (30) is lower than 

2 el^°lal^'^'l2l^"l. 
Ge?ii(y(+)) 

And so, performing the same calculations as for the proof of Lemma 2, we obtain, for 
the second part of (28), 



ILL 

xeZd y^Px y(+)Cy\{yo} 
<Cinv||<p|icT", 



ri(lF(yW,or*'''''-'')^£(r(-'.o)^ror^5<- 



(31) 



where 



1 - [2'?(i?2 +2i?)]2£ 1 - [2'?(i?2 +27;)]-2e ^^^v 
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Inserting (29) and (31) in (28), we conclude: 

rX(p)-M,^((p)|<c||(p||cT", 

where we renamed C + Qnv to C for simplicity. 



7 Exponential Decay of Correlations 

We will now prove a well-known consequence of Theorem 1 : the invariant measure 
A'i'nv presents exponential decay of correlations in space. 

Proof (Corollary) Since 5'^' belongs to ^'+'(1,0), Theorem 1 implies that for some 
C <oa and some d < 1, we have, for all ^ G N, 

r*5W((pv/)-<'!(<Pr)| <C||(pv/|lcT*; 
Tk5^-\<p)-^^{<p)\<C\W\\ah 

But 5'^' is a product measure and the interactions are local, so 

r*^5<+'(9r) = r'^5'+'((p)r'^5'+'(v/) 

as long as 2A:max„gj/ < d{(p,\i/), where ||-||[ is the Manhattan norm. Since 
d{(p,\f/) > 0, we also have 

llwll < !I<pIIIv^L + |(pL||v^|| <°o. 

Hence, as long as 2A-max„gy < d{(p, i//), 

-Mi^!(<p)Mi^:(r) < 2C{M \wL + \(pL \M)<y'- 



If wetakefc= \d{(p, v/)/2max„gf/ ||m||j] - 1 and choose C' = ^ and 77 = , 
we obtain the desired inequality. 

The invariant measure ^^'^l also exhibits exponential decay of correlations in time. 
Proof (Theorem 2) Theorem 1 applied to 5'^' € ^'^'(1,0) implies that, for all A: G N, 

^il2i<pT»-T'8^^\(pT"xi,)\<C\\(pT"xtf\\a'. (32) 

Now, remembering definition (2) of the semi-norm ||.||, we notice that 

||(prv/||<2|(pUiv/U|supp((pr»|. 

But the conical structure of the space-time influence of spins implies 



\supp{(pT" < |supp((p)| + (diam(supp(v/)) +2nmaxj|«||j )'' 
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SO there exists a finite constant ^ such that ||^r"\//|| < ^«'' for all n G Nq. 
Theorem 1 also yields 

Mi^('pR"v('/)-Mi^v('P)7^'+"5"(v^)| <c|<PLIlrl|(T'^+"; (33) 

Mi^i(9)r*^+''5W(v/)-r'^5W((?))r*+"5'+'(v/) <C|v/L||(P||(t^ (34) 
In order to find an upper bound for ( ((p - r^'5'+'(9)) T"\^)\, we rewrite it 



as 



(y/) , where the signed measure /i^*' e (X) is defined by 



/i^ is a measure of zero mass so it satisfies (15): 



n 1 



(=1 



The last expression is similar to the second term of (28) with /x^*^' instead of 

and we can then perform the same argument as in the proof of Theorem 1 in order to 



bound 



extending again the collections of Toom graphs up to time —k. 

The only change in these calculations is the presence of an extra factor <p — (^) 
whose supremum norm is bounded by 2 |^|^. Provided we keep track of this factor, 
the calculations which lead to (31) still hold here: 



(35) 



Combining (32), (33), (34) and (35) and taking the k ^ °o limit ends the proof. 



Theorem 2 implies that ji^^^^ is not only ergodic but also mixing, that is, for any 



continuous functions (p, y/, we have 



lim 



<":('pr"v/)-<^:((p)M«(v/) 



= 0. 



Indeed, by definition of^{X), the set of continuous functions with finite support is 
dense m^{X). 
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